Close encounters of three black holes 
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We present the first fully relativistic longterm numerical evolutions of three equal-mass black holes 
in a system consisting of a third black hole in a close orbit about a black-hole binary. We find that 
these close-three-black-hole systems have very different merger dynamics from black-hole binaries. In 
particular, we see complex trajectories, a redistribution of energy that can impart substantial kicks 
to one of the holes, distinctive waveforms, and suppression of the emitted gravitational radiation. 
We evolve two such configurations and find very different behaviors. In one configuration the binary 
is quickly disrupted and the individual holes follow complicated trajectories and merge with the 
third hole in rapid succession, while in the other, the binary completes a half-orbit before the initial 
merger of one of the members with the third black hole, and the resulting two-black-hole system 
forms a highly elliptical, well separated binary that shows no significant inspiral for (at least) the 
first t ~ lOOOM of evolution. 

PACS numbers: 04.25.Dm, 04.25. Nx, 04.30.Db, 04.70. Bw 



Introduction: The recent dramatic breakthroughs in 
the numerical techniques to evolve black-hole-binary 
spacetimes P, 0, [3| has led to rapid advancements 
in our understanding of black-hole physics. Notable 
among these advancements are developments in mathe- 
matical relativity, including systems of PDEs and gauge 
choices '3, 'H', 'i'], the exploration of the cosmic censor- 
ship [7,, ^, J,, dlKand the application of isolated hori- 
zon formulae @, la 111, [11 [11 tii. These breakthroughs 
have also influenced the development of data analy- 
sis techniques through the matching of post-Newtonian 
to fully-numerical waveforms (igI . 
the recent discovery of very large 
0, H [m, [^, 23, ' 24i,|2i,|2r 
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Similarly, 
merger recoil kicks 
has had a great impact 
in the astrophysical community, with several groups now 
seeking for observational traces of such hig h s p eed holes 
as the byproduct of galaxy coUisions [271, [2^. In this 
letter, we continue our quest to discover new astrophysi- 
cal consequences of black- hole interactions by simulating 
close encounters of three black holes to see the differ- 
ent behaviors introduced by the finite size of the holes, 
their nonlinear interactions, and the radiation of gravi- 
tational waves, as described by General Relativity. We 
find that the three-body relativistic problem shows far 
richer dynamics than the two-body problem, akin to the 
rich three-body dynamics in Newtonian gravity, but with 
added complexity due to mergers. 

Three-body and four-body interactions are expected 
to be common in globular clusters 2^, S^l , and in galac- 
tic cores hosting supermassive black holes (when stellar- 
mass-black-hole-binary systems interact with the Super- 
massive black hole) . Hierarchical triplets of massive black 
holes might also be formed in galactic nuclei undergo- 
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ing sequential mergers [31|, |32] ■ The gravitational wave 
emission from such systems was recently estimated using 



post-Newtonian techniques [33| . 

Techniques: We evolve the three-black-hole data-sets 
using the LazEv [34| implementation of the 'moving 
puncture approach' [1[3|. We use the Carpet [35l| driver 
to provide a 'moving boxes' style mesh refinement. In this 
approach refined grids of fixed size are arranged about the 
coordinate centers of each hole. We use AHFinderDi- 
RECT js^l to locate apparent horizons. We extract the 
waveform on spheres centered about the origin and ex- 
trapolate the radiated energy/momentum to r = oo (the 
waveforms do no change qualitatively for r > 50±20Af). 

Results: We chose one configuration (3BH1) with 
purely ad-hoc momentum parameters, which merged rel- 
atively quickly, to test the convergence and accuracy of 
our code. The initial data parameters for these config- 
urations are summarized in Table H] We evolved these 
configuration using 11 levels of refinement and a finest 
resolution of ft. = M/80. The outer boundaries were lo- 
cated at 640M. In addition we evolved the 3BH1 config- 
uration with grid-spacings rescaled by 5/6 and (5/6)^ to 
test convergence. In the table, the horizon mass is the 
Christodoulou mass, where — \jTn^„ + <S'2/(4m? j.), 
S is the magnitude of the spin of the hole, and m\„ is 
the irreducible mass. 

It is interesting to note that the same techniques used 
for black-hole binary evolutions work for configurations 
of three (and, according to a brief test by the authors, at 
least 22) black holes. We tested the convergence of our al- 
gorithm with three black holes by evolving configuration 
3BH1 with three resolutions (Af/80, Af/96, Af/115.2). 
We chose this configuration since it merges relatively 
quickly, thus reducing the computational expense. The 
resulting waveform converges to fourth-order. Note that 
late-time accumulation of errors can have a significant 
effect on the trajectories of 3-black-holc systems due to 
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TABLE I; Initial data parameters. {xi,yi,0) and (p^,pj',0) 
are the initial position and momentum of the puncture i, 
is the puncture mass parameter, and mf is the horizon mass. 



Config 


3BH1 


3BH101 


3BH102 


xi/M 




-3 52462 


-3 52238 


7/1 /M 


2.23413 


2 58509 


2 58509 


nf /M 
v"JM 
m\IM 


-0.0460284 


-0.0782693 


0.0782693 


-0.0126181 


-0.0400799 


-0.0433529 


0.315269 


0.318143 


0.317578 


mf /M 


335555 


336201 


335721 


X2/M 


-2 40856 


-3 52462 


-3 52462 




-2.10534 


-2 58509 


-2 58509 




0.130726 


0.0782693 


-0.0782693 


pI/m 

ml/M 


-0.0126181 


-0.0400799 


-0.0433529 


0.315269 


0.318143 


0.317578 


/M 


0.3405205 


0.336241 


0.335767 


Xi/M 


4.8735 


7.04923 


7.04476 


Vi/M 


0.0643941 








Pl/M 


-0.0846974 








pI/m 

ml/M 


0.0252361 


0.0801597 


0.0867057 


0.315269 


0.320815 


0.318585 


m^/M 


0.332198 


0.333115 


0.331270 



5 - 







-5 




_10 I ' ' • ' ' 1 

-20 -15 -10 -5 5 10 

x/M 

FIG. 1: The horizon trajectories for configuration 3BH101. 
The three black holes are initially located at the points la- 
beled by BHl, BH2, and BH3, respectively. BHl and BH2 
form a quasi-circular binary, which is disrupted by BH3. BH2 
and BH3 merge at point Ml. The BHl and the BH2-BH3 
merger remnant continue to orbit each other throughout the 
simulation. 
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their inherent sensitivity to changes in configuration. We 
have confirmed, through the use of a 8th-order accurate 
code, that the trajectories presented here show the cor- 
rect quahtative behavior. 

We now show results for two similar initial configu- 
rations with qualitatively different outcomes. To aid in 
the discussion, we will denote the two holes in the bi- 
nary with BHl and BH2, and the third hole with BH3, 
where BHl is initially located at y > 0. We determined 
our initial data parameters by choosing a fiducial binary 
configuration with orbital frequency Mb^b = 0.04 and 
angular momentum Jb/M^ — 0.9104975. We then treat 
the binary as a point particle of mass Mb = 2/3M and 
spin ub/Mb = 0.9105 along with a non-spinning point 
particle of mass M/3, and choose position and momenta 
parameters such that this two particle system is in a 
quasi-circular orbit (up to 3 PN) at a separation equal to 
twice the binary's separation. We set up the systems so 
that the binary's orbital angular momentum is aligned 
with the total orbital angular momentum (configuration 
3BH101), and anti-aligned (configuration 3BH1G2). Con- 
figuration 3BH101 has BH2 and BH3 merging after the 
binary completes nearly a half of an orbit. The result of 
this interaction is to significantly push BHl away from 
the merger remnant, producing a new, highly elliptical, 
binary, with large orbital separation oscillating in time 
from 10.4M to 23.5M (see Fig. H]). The 3BH101 wave- 
form (Fig. [2]) shows a burst of radiation from the BH2- 
BH3 merger, as well as a small pulse at t ^ 700 Af which 
corresponds to the point of closest approach of the BH2- 
BH3 merger remnant with BHl. We stopped the evo- 
lution a,t t ^ lOOOM due to computational expense and 
boundary contamination. 
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FIG. 2: The (Z = 2,m = 2) mode of i>4 for 3BH101. The 
BH2-BH3 merger waveform is centered at f ~ 185M, the 
small pulse at f ~ 700 was produced by the close approach of 
the BH2-BH3 merger product to BHl. 



Configuration 3BH102 displays very different behav- 
ior, as seen in Fig.[3l Here the binary is disrupted almost 
immediately, and the individual holes follow complicated 
trajectories (note that the trajectories are similar to the 
Greek letters 7, r, and a). BH3 and BHl merge when 
BH3 almost completes 1.25 orbits. The BH3-BH1 merger 
product then quickly merges with BH2. The resulting 
waveform shows a double merger as seen in Fig. |4l It 
is important to note that the dramatic difference in dy- 
namics between 3BH101 and 3BH102 is not strictly due 
to corotation versus anti-corotation of the binary, but 
rather is a result the precise configuration of the binary 
when the third black hole approaches. 

For all three configurations the radiated energy and an- 
gular momenta were a fraction of that for a quasi-circular 
equal-mass binary. This is due to the grazing type merg- 
ers, as also seen by the suppression of the radiated an- 
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TABLE II: Total radiated energy, momentum, and angular 
momentum, as well as horizon mass spin and merger time, 
for 3BH102 and 3BH1. (* 3BH101 did not merge in the time 
alloted and we only report values for the first merger). 



Config 



-2.5 2.5 
x/M 



£;,ad/M(xlOOO) 
J^lJM^x WO) 
P^l^/M{kms-^) 

Mu/M 

tmi/M 

tM2/M 



3BH1 



3BH10r 



6.06 ± 0.02 
2.92 ±0.01 
20.0 ± 1.9 
-22.9 ±2.4 
0.9835 

0.532 

~ 27 

~ 40 



4.4 ±0.8 
3.6 ±0.5 
50 ±40 
27± 13 

+ + + 
* + + 

~ 115 



3BH102 



6.1 ±0.7 
4.4 ±0.2 

-1.3 ± 14 
-15± 13 
0.9885 
0.465 
~ 218 
~ 280 



FIG. 3: The horizon trajectories for configuration 3BH102 
with ticks every 45M of evolution. The three black holes 
are initially located at the points labeled BHl, BH2, and 
BH3, respectively. BHl and BH2 form a quasi-circular bi- 
nary, which is almost immediately disrupted by BH3. BHl 
and BH3 merge at point Ml, and then merge with BH2 at 
M2. 
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FIG. 4: The = 2, m = 2) mode of ^4 for 3BH102 shows two 
merger waveform signals arising from merger Ml and M2. 



gular momentum. We expect that the fuU waveform for 
3BH101 wiW show significantly more radiation as the sys- 
tem eventually 'circularizes' and merges. 

We calculate the mass and spin of the first merger 
remnant and the final remnant using the fitted exponen- 
tial decay rate and frequency of the quasi-normal behav- 
ior (33. For 3BH102, we fit the real and imaginary parts 
of the {£ — 2,m = 2) mode separately and find an /Mh = 
0.479,0.378 and Mh = 0.716,0.678 from fits of the real 
and imaginary parts of "04 for the first merger product, 
and qh/Mh = 0.478,0.580 and Mh = 0.994,1.04 from 
fits of the real and imaginary parts of V'4 for the final 
merger remnant. Note that the imaginary part of "04 pro- 
vides a better estimate (i.e. closer to the expected 0.66) 
mass for the first remnant, but the real part provides a 
better estimate of both a and Mh for the final remnant. 
We note the quasi-normal frequencies were lo = 0.64, 0.46 
for the first and final merger remnants respectively. For 



3BH101 we were only able to fit the real part of -04 due to 
significant contamination from other modes in the imag- 
inary part. We find a narrow region of width lOM where 
the waveform shows nearly exponential decay. A fit to 
this region yields u = 0.63, with a corresponding mass 
and spin of Mh = 0.665 and a/Mn = 0.293. The rem- 
nant masses, spins, and merger times are given in Table HIl 
for the 3BH1 and 3BH102 configurations. The numbers 
quoted above should be taken only as indicative of the 
expected values of these parameters. Further runs, at 
higher resolutions, will be needed in order to establish 
the errors in these values. 

Discussion: The relativistic study of quasi-circular or- 
bits of a binary in the presence of a third comparable- 
mass hole, as an initial value problem, was studied in [sst . 
Here, by studying their dynamical evolution we find that 
the third hole perturbs the system to the extent that 
no true binary orbit is seen. We found that the close 
encounter of this third body can both trigger a quick 
merger of the three-body system, as well as impart a sig- 
nificant kick to one of the holes, producing a new long- 
lived, highly-elliptical binary. The generic effect of the 
third black hole is to reduce the gravitational radiation. 
This happens for two reasons. First, close-three-body in- 
teractions lead to grazing collisions, which emit far less 
radiation than quasi-circular mergers. Second, the result- 
ing binary orbit will be elliptical, which is less efficient 
at emitting gravitational radiation than circular orbits at 
the final stages. Note however, that although we report 
radiated energies that are 1/5"^ that for a typical binary, 
here we scale the energy by the total mass. If we scale 
the radiated energy with the initial binary's mass, then 
the rescaled radiated energy would be 3/2 times larger. 
The close-three-body systems also appear to be shorter 
lived than typical binaries. 

The three-black-hole waveforms (See Figs. [2] and [4]) are 
distinct from the robust and simple form of the binary- 
black-hole waveform © 113, EH, El. In addition, there 
seems to be a large exchange of energy among the com- 
ponents of the triple system, which occurs on a much 
shorter timescale than the radiation. It is important to 
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note that these three-black-hole interactions provide a 
mechanism for producing highly-elliptical close-binaries, 
which would otherwise have circularized (due to emission 
of gravitational radiation during the inspiral) . 

Investigations of the Newtonian encounters of three 
bodies show that such encounters generically lead to the 
breakup of the system into a binary and the third body 
that escapes [43| in a 'water-shed effect'. The distribu- 
tion of the eccentricity of the remaining binar y is bell 
shaped around e = 0.3 for compact systems |43 . |45| . 
Classical studies Q show that the probability of ex- 
change of the binary companion in a triple system is sur- 
prisingly high for all comparable masses, reaching near 
one for more massive ma [47| . The motion of the system 
can be chaotic, due to small denominators. The finite 
size of the black holes represents a natural regularization 
to the problem, and the dissipative effects of the gravi- 
tational radiation can prevent some configurations from 
becoming chaotic. 

We have found that 3-black-hole systems exhibit com- 
plicated orbital dynamics analogous to the rich 3-body 
Newtonian dynamics, but with the added complexity in- 
troduced by mergers. Further study of this problem, in- 
cluding a comparison of Newtonian and relativistic dy- 
namics, will be reported in a forthcoming paper, where 
we also examine configurations where the triple system 
is disrupted. 
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